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From proper motions to star cluster dynamics: measuring 
velocity dispersion in deconvolved distribution functions 
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ABSTRACT 

We investigate the effect that the usuaUy large errors associated with ground-based 
proper motion (PM) components have on the determination of a star cluster's velocity 
dispersion (ay). Rather than histograms, we work with PM distribution functions 
(PMDFs), taking the la uncertainties formally into account. In this context, a cluster's 
intrinsic PMDF is broadened by the error distribution function (eDF) that, given the 
average error amplitude, has a width usually comparable to the cluster PMDF. Thus, 
we apply a Richardson-Lucy (RL) deconvolution to the PMDFs of a set of relatively 
nearby and populous open clusters (OCs), using the eDFs as point spread functions 
(PSFs). The OCs are NGC1039 (M34), NGC2477, NGC2516, NGC2682 (M67), 
and NGC7762. The deconvolved PMDFs are approximately Gaussian in shape, with 
dispersions lower than the observed ones by a factor of 4-10. NGC 1039 and NGC 2516, 
the nearest OCs of the sample, have deconvolved ay compatible with those of bound 
OCs of mass - IO^Mq. NGC 2477 and NGC 2682 have deconvolved PMDFs with a 
secondary bump, shifted towards higher average velocities, which may be an artefact 
of the RL deconvolution when applied to asymmetric profiles. Alternatively, it may 
originate from cluster merger, large-scale mass segregation or, least probably, binaries. 

Key words: ( Galaxy:) open clusters and associations: general 



> 1 INTRODUCTION 



1 



Star clusters continually undergo mass segregation and 
evaporation, tidal interactions with Galactic substructures, 
shocks with giant molecular clouds, as well as mass loss due 
to stellar evolution. By decreasing the total cluster mass 
- and the collective gravitational potential, these processes 
affect the internal dynamics and accelerate the cluster dy- 
namical evolution. Thus, the escape velocity (vesc) and, to a 
lesser degree, the space velocity dispersion (o"i,), are expected 
to continually change with time. As a consequence, the ma- 
jority of the open clusters (OCs) dissolve in the Galactic 
stellar field long before reaching an age of ^ 1 Gyr (e.g. 
iLamers et al]|2005l : ICoodwin fc Bastianl 120061 ). In this con- 
text, the dynamical evolution of a star cluster depends crit- 
ically on the balance between ct„ and Vesc- 

For a cluster that is approximately in virial equilib- 
rium, the space velocity dispersion can be computed as 
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l|Spitzeilll987l '). where G is the gravitational constant, Md 
is the dynamical mass (assumed to be stored only in sin- 
gle stars), 77 ~ 9.75 is a constant, and Reff is the effective, 
or projected half-light radius. As a scaling factor, bound 



clusters with Md ~ 10 M0 and Reff ~ Ipc are expected 
to have (t„ ~ lkms~^. Since the majority of the Galactic 
OCs have masses lower than 10^ Mq (e.g. iPiskunov et al.l 
[2003), velocity dispersions lower than, or of the order of 
av ~ lkms~^ should be a rule. 

Obviously, the above expectation does not apply to clus- 
ters that are on their way to dissolution, either at the earliest 
evolutionary stages (less than a few 10 ^ yr) or much later, a t 
the remnant phase (several 10* yr, e.g. IPavani fc Bicall2007i '). 
Such a super-virial state in very young clusters is related 
to the impulsive expulsion of the residual parental molecu- 
lar cloud gas, driven prim arily by the strong winds o f mas- 
sive stars and supernovae ICoodwin fc Bastianl lioo^ l. thus 
leading to the high dissolution rate of young clusters (e.g. 
iLada fc Ladall2003l l. Examples of dissolving young OCs with 
velocity d ispersions that by far exceed the expected dynami- 
cal value (|de Griis et al.|[2008i ) are NGC 2244, with an age of 
^ 5Myr and ^ 35k ms~^ (|Chen. de Griis fc Zhao|[20oi : 
iBonatto fc BicalbOOQbl l. Cr 197 an d vdB92, both with an 
age of ~ 5 Myr and cr^ ~ 20 km s"^ (|Bonatto fc Bicall2010al ') 

A first-order estimate of a star cluster's (t„ can be 
obtained by measuring the proper-motion (PM) distribu- 
tion of the member stars. Having obtained the PM com- 
ponents in right ascension (/Iq cos((5)) and declination (us) 
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of each member, the aa = cos((5) ^^'^ '^^ ^ "^A** 
persions can be estimated from the shape of the respec- 
tive distributions. Ahernatively, if the systemic PM com- 
ponents are known, one can compute directly the tangen- 
tial or on-sky PM po, 
sion a OS = '/o. 



75 



take cr„^ = |(cr, 



= \/ Ha cos(5)^ + fj.s'^ and disper- 
Then, assuming isotropy, we can 
At this point, it should 



be mentioned that, contrary to radial velocity measure- 
ments, proper motions are essent ially unaffected by binarity 
l|Kouwenhoven fc de Gr"ii3 120081 ). This occurs because ra- 
dial velocities are instantaneous measurements, while proper 
motions involve a (usuaUy larg e) timespan. In this sense, 
iKouwenhoven fc de GriisI (|2008l ) show that for clusters of 
mass ~ 1000 Mq , binaries produce an observed radial veloc- 
ity dispersion several times higher than the dynamical one, 
thus leading to an artificially-overestimated cluster mass. 

Comprehensive databases such as the Third U.S. 
Naval Observa t ory C CD Astrograph Catalog (UCAC3, 
IZacharias et al.l |2010| ) and the The Naval Observa tory 
Merged Astrometric Dataset (NOMAD, Zacharias c t al.l 
(2003) provide PM components for huge amounts of stars, 
either isolated or in clusters. However, a drawback of such 
all-sky, ground-based PM databases is the fact that the la 
uncertainties in jia cos((5) and jis are usually large. The effect 
that such errors have on the determination of (t„ , which usu- 
ally involves building PM histograms with bin sizes smaller 
than the average la errors, is a constantly neglected point. 
Instead of histograms, our approach in this paper is to work 
with PM distribution functions (PMDFs), in which the er- 
rors are formally taken into account. 

This paper is organised as follows. In Sect. [2] we build 
PMDFs and compare them with the classical histograms. 
In Sect. [3] we select some relatively nearby and populous 
OCs as test cases. In Sect. |4] we discuss the deconvolution 
method that we apply on the PMDFs. In Sect. [5] we present 
the results obtained so far and discuss some observational 
limitations. Concluding remarks are given in Sect. (6] 



2 PROPER MOTION ANALYSIS: 

HISTOGRAMS OR DISTRIBUTION 
FUNCTIONS? 

When dealing with PMs, the first step is to isolate the (most 
probable) cluster members, subtract the systemic compo- 
nents of /ic cos((5) and from each star, and then compute 
the tangential or on-sky PM fj,oB. The systemic components 
can be directly computed from the data or obtained in OC 
databases (e.g. WEBDAQ or DAMLofl). 

Usually, the next step involves building PM histograms 
with a given bin width to estimate the dispersion either of 
the components or the on-sky. However, except for a few par- 
ticular cases with high-quality astrometrj|j, ground-based 
PMs usually have significant errors, in general larger than 



^ www.univie.ac.at/webda 

Catalog of Optically Visible Open Clusters and Candidates 
www.astro.iag.usp.br/ wilton 

^ For instance, the globular cluster NG C 6397 with Hubble Space 
Telescope WFPC2 data from 2 epochs llRicher et al.ll200Sl) . 
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Figure 1. Histogram and PMDF (heavy-solid line) of NGC 2682 
(top panel) and NGC 2477 (bottom) for the stars located within 
-R < 20'. The amplitude of the maximum, average, and minimum 
Icr errors are illustrated (horizontal lines). Also shown is the error 
distribution function (dashed line). 



the adopted histogram bins. Thus, by ignoring the error am- 
plitude, histograms should be taken only as first-order rep- 
resentatives of a cluster's intrinsic PM distribution. Wider 
bins would obviously minimise this effect, but would also 
degrade the PM resolution. The net result would be arti- 
ficially high values of the dispersion. This is illustrated in 
Fig. [T] for the cos(5) component of the OCs NGC 2682 
and NGC 2477 (Sect. [3}, where we show the classical his- 
togram built with bins of size Afia cos{S) = 1 mas yr~^. Also 
shown are the average, minimum, and maximum values of 
the la errors. Clearly, in both cases the average la error 
corresponds to approximately the core (« 4 bins) of the his- 
togram. It should be noted as well the introduction of a sort 
of high-frequency noise, which originates from the Poisson 
fluctuation associated with the bin size. 

Thus, instead of histograms, the correct way of dealing 
with the above issue is by explicitly incorporating the PM 
component errors in PMDFs, which are defined as the frac- 
tional number of stars per interval of PM, PMDF = <j>{fi) = 
dN /d^, where ^ is any component among /i^ cos(5), fis, and 
yLoa ■ The errors are incorporated in the PMDFs by assuming 
that they are normally (i.e. Gaussian) distributed. Accord- 
ingly, if measurements of a given parameter x are normally 
distributed around the average x with a standard deviation 
e, the probability of finding it at a specific value x is given 



by P(x) 



i(^) 



Wc start by defining a set of PM bins spanning the 
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whole range of values observed in /iQCOs((5), fis, and fios- 
The bins have smaller widths around the PMDF peak than 
for higher or lower velocities (for preserving profile resolu- 
tion - App. [XJ. Then, for a star with PM components and 
la uncertainties fia cos{5) ± e^^ cos(S) ' ^ ''M* ' A^os ± <^Mos > 
we compute the probability of the PM of that star to be 
in a given bin, which is simply the difference of the error 
functions computed at the bin borders. By doing this for all 
stars, we end up with the number-density of stars in each PM 
bin, the integral of which over the whole range of PM val- 
ues is simply the number of stars. The Ha cos{5) PMDFs of 
NGC 2682 and NGC 2477 are shown in Fig. [T] By construc- 
tion, the PMDFs are much smoother than the histograms 
(the high-frequency noise has been naturally removed) and, 
because of the broadening effect due to the error spreading 
procedure, they are also somewhat wider and have a lower 
amplitude. Both the histograms and PMDFs are definitely 
non-Gaussian, especially because of the broad wings. 

Now, the same procedure is applied to the uncertain- 
ties, thus resulting in the intrinsic error distribution func- 
tion (eDF). As anticipated by the amplitude of the aver- 
age la error, the eDF has a width comparable to that of 
the PMDF for both NGC 2682 and NGC 2477. In this con- 
text, the eDF plays a role of a PM point-spread function 
(PSF), which tends to broaden the intrinsic cluster PMDF 
on a degree that depends essentially on the PMDF and eDF 
widths. Thus, our approach here is to deconvolve the ob- 
served PMDF of selected OCs (Sect. [3]) using the intrinsic 
eDF as PSF (Sect.lJ). 



3 TEST CASES 

As test cases we searched for OCs that are located rela- 
tively nearby (for allowing detection of low-PM components 
- Sect.[Sl), away from central directions and the disk (to min- 
imise field-star contamination), and with a wide range of 
ages. Additionally, the candidates should have a significant 
number of stars with available PM components in UCACvQ 
(for more representative PMDFs). The sample of OCs meet- 
ing our criteria are NGC 1039 (M 34), NGC 2477, NGC 2516, 
NGC 2682 (M 67), and NGC 7762; their fundamental param- 
eters are listed in Table [1] In short, the selected OCs have 
ages within 0.2 - 4 Gyr and distances from the Sun within 

~ 0.4 - 1-3 kpo 

2MAS#I (|Skrutskie et aLlbOOd ) colour-magnitude dia- 
grams (CMDs) extracted from the central region (R < 5') 
are shown in Fig.O The 5' boundary was selected for practi- 
cal reasons. Our analysis depends on the number of member 
stars in two different cluster regions (Sect.|4]), and the stel- 
lar density profiles (Fig. ^ show that R — 5' represents a 
compromise between adequate statistics and a change from 
RDPs following a pure power-law {R > 5') to a power- 
law fiattened by a core (R < 5'). For NGC 26 82 (M67) we 
also use the ground-based CCD astrometry of lYadav et al.l 



* Main features of UCAC3 are: complete sky coverage, merging 
of several PM catalogues, new data reduction with reduced er- 
rors, significantly improved photometry from CCD data, etc. See 
I http:/ / cds arc.u-strasbg.fr/ viz- bin/ Cat ?!/ 315\ for further details. 
" The Two Micron All Sky Survey, All Sky data release - 
\http://www.ipac.caltech.edu/2mass/releases/allsky/\ 
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Figure 2. Extracted within R < 5', the CMDs have been built 
with 2MASS photometry, exce pt for pane l (b) th at was built with 
the VI CCD photometry of Y adav et al.l ll2008h . Only stars that 
occur within the colour-magnitude filter (shaded polygon) are 
considered in the analyses. 



(|2008| ) as an independent data set for cr oss-checking our re - 
sults. The PM components provided bv lYadav et al.l (|2008l l 
are already corrected for the cluster's systemic velocity. The 
2MASS photometry was extracted from ViziePlfl in a wide 
circular field of radius i?ext = 90', which is adequate for 
determining the cluster extension and background level (see 
below). We also build a colour-magnitude filter for each clus- 
ter, which is used to iso late probable member stars (see, e.g. 
iBonatto fc Bicall2010bl for a discussion on the use of such 
filters). Only stars with colour and magnitude compatible 
with the filter are used in the subsequent analyses. Age, 
distance from the Sun, and reddening are derived by fit- 
ting solar-metallicity Padova is ochrones (jCirardi et al H2OO2I : 
iBonatto. Bica fc GirardU 120041 ') to the CMDs. The derived 
fundamental parameters (Table [T| are similar to those in 
WEBDA. 

After isolating the probable member stars with the 
colour-magnitude filters, we use them to build the stellar 
radial density profile (RDP), which is the projected stellar 
number density profile around the cluster centre (Fig. |3]). 
Note that the cluster centres have been computed by an al- 
gorithm that searches for the highest stellar density in the 
innermost bin and, at t he same time, the smoothest RDP 
IIBonatto fc Bical l2010al ). Working with colour-magnitude 

\http://vizier.u-strasbg.fr/viz-bin/VizieR ?-source= II/246\ 
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Table 1. Fundamental parameters and cluster systemic PM components derived in this work 



Literature This work 



Cluster 


a(2000) 


(5(2000) 


£ 


b 


Age 


Av 


dQ 


Rrup 


5' 


(/^Q, cos((5))sj/3t 






(hms) 


{° 


"') 


n 


n 


(Myr) 


(mag) 


(kpc) 


(pc) 


(pc) 


(mas yr"-') 


(mas yr~-^) 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


NGC 1039 


02:42:45 


+42 


45 


42 


143.66 


-15.61 


300 ± 50 


0.2 ± 0.2 


0.49 + 0.03 


4.3 + 0.7 


0.7 


-0.55 + 0.01 


-5.81 + 0.01 


NGC 2477 


07:52:10 


-38 


31 


48 


253.56 


-5.84 


900 ± 100 


0.8 + 0.2 


1.31 + 0.07 


13.3 + 1.9 


1.9 


-0.83 + 0.07 


+1.89 + 0.06 


NGC 2516 


07:58:04 


-60 


45 


12 


273.82 


-15.86 


200 ± 50 


0.3 ± 0.2 


0.38 + 0.04 


3.3 + 0.3 


0.6 


-5.60 + 0.07 


+10.74 + 0.03 


NGC 2682 


08:51:18 


+ 11 


48 


00 


215.70 


+31.90 


4000 ± 300 


0.0 + 0.1 


0.79 + 0.05 


6.9 + 0.7 


1.2 


-8.61 + 0.02 


-4.92 + 0.05 


NGC 7762 


23:50:01 


+68 


08 


18 


117.22 


+5.85 


2000 ± 500 


1.4 + 0.4 


0.96 + 0.14 


7.0 + 1.4 


1.4 


-4.67 + 0.01 


-2.53 + 0.03 



Table Notes. Cols. 4 and 5: Galactic coordinates; Col. 8: distance from the Sun; Col. 9: cluster truncation radius (in pc); Col. 10: R = 5' 
in absolute units; Cols. 11 and 12: deconvolved systemic proper motion components, computed with the stars within R < 20'. 
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Figure 3. Stellar RDPs of the selected DCs. The central (0 < 
R(') < 5), outer (5 < i?(') < 20 or 5 < R{') < 15 - top right), 
and field (F) regions are indicated. The field contamination level 
towards the central parts is illustrated by the shaded polygon. 



filtered photometry minimises field contamination and en- 
hances the RDP contrast with the fore/background. The 
position (and error) of each RDP point along the R axes in 
Fig. |3] corresponds to the average (and la uncertainty) dis- 
tance to the cluster centre of the stars within each bin. We 
also estimate the cluster truncation radius (-Rrdp), which 
is simply the distance from the cluster centre where RDP 
and background are statistically indistinguishable (Fig. [3]). 
On average, the boundary R = 5' represents Ri 1/6 of -Rrdp 
(Table m. The near- infrared RDPs show that the selected 
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Figure 4. RL deconvolution applied to the fia cos((5) (left panel), 
lis (centre), and /los (right) PM components for the stars located 
in the region R < 20' of NGC 2682. Top panels show the observed 
cluster (solid line) and field (shaded area) PMDFs, and the PSF 
(dotted line). Bottom: field-subtracted (shaded area) and decon- 
volved (solid line) PMDFs. 



P Cs have an extensi on reaching ~ 20— 30' (the optical data 
of lYadav et al ] |2008l are restricted to ~ 15'). They also show 
that field stars with the same colour and magnitude as the 
probable members are still present in the central region, in 
varying amounts for the different OCs. The residual field 
contamination will be taken into account in the PM analy- 
sis (Sect. 131). 



Proper motion components were obtained in UCAC3 
(that also provides the 2MASS photometry for each star) 
based on the same central coordinates and extraction radius 
as those used for building the CMDs and RDPs. Afterwards, 
we applied the respective colour-magnitude filters (Fig. [2]) 
before computing the PM distributions (see below). 



4 PMDF DECONVOLUTION 

By construction, the observed PMDFs are broadened by a 
PSF that corresponds to the eDF (Sect.[2|. Thus, to uncover 
the intrinsic PMDFs, we apply the iterative Richardson- 
Lucy (RL) decon volution method proposed by iRichardsonI 
l|l972l ') and ILucvI (1974). The RL deconvolution conserves 
the PMDF integral (in the present context the number of 
stars), but has a relatively slow convergence rate. 

We illustrate the RL deconvolution with the cos (5), 
lis, and iJios components of the OC NGC 2682 (Fig.[4|. As the 
first step we define the PM-bin size distribution fApp. [A]), 
and compute and subtract the systemic PM components 
(Table [l| from ^aCos{5) and ^s. Actually, this is an iter- 
ative step, in which we first subtract the field contamina- 
tion, deconvolve the resulting PMDF and obtain the sys- 
temic PM. We use the region R < 20' in the analysis as a 
compromise between the number of member stars and con- 
trast with the field level (Fig. [3|. Then we build the cluster 
and field PMDFs, together with the eDF, which has a width 
similar to that of the cluster PMDF. Next, the contaminant 
PMDF of the residual field stars, which is especially seen in 
fjLaC0s{5) as a bump centred at ~ 6 mas yr~^, is subtracted 
from the observed PMDFs. Finally, the RL deconvolution is 
applied to the field-subtracted PMDF. 

Besides the slow convergence rate, the RL deconvolu- 
tion is also known for not having a universal convergence 
criterion. Thus, for determining the number of deconvolu- 
tion iterations (Nu), we simply sum the squared difference 
between successive iterations over all bins of the deconvolved 
PMDF, defined as 

X? = J2 {PMDF{xi)j - PMDF(x,),.^_Y, 

i 

where j is the current iteration and Xi is the i"* bin along 
the PM axis (either /iaC0s(5), /if, or /ios), and compute 
the fractional variation of x, /x — 1 ~ Xj/Xj-i- After a 
series of tests to check changes in deconvolution parame- 
ters with the number of iterations, we arbitrarily decided to 
stop wheniVit = 200 (or = 0.003-0.005; App.^J. Under 
this criterion, the deconvolved PMDFs end up with signifi- 
cantly lower dispersions than the observed ones (Figs. [4] to 
[6] and Table [2]), without significant added noise or artifacts 
(App. Besides, contrary to the observed PMDFs, the 
shape of the deconvolved ones is approximately Gaussian. 



5 DISCUSSION 

Similarly to NGC 2682 (Sect.[4|, the RL deconvolution was 
applied to the field-subtracted PMDFs of the remaining 
OCs. Here we present the results of this procedure but, to 
avoid redundancy, we restrict the discussion to the on-sky 
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NGC2682 (UCAC3) NGC2682 (optical) NGC1039 




(mas yr" ) 

Figure 5. Field-subtracted (shaded area) and deconvolved (cir- 
cles) on-sky PMDFs built for different cluster regions: central 
(top), outer (middle) and overall (bottom). The solid line shows 
the Gaussian fit to the deconvolved PMDF. The number of mem- 
ber sta rs in each region is i ndicated. The PMDFs with the optical 
data of lYadav et~aU l|2008h for NGC 2682 (central panels) are re- 
stricted to _R < 15'. 



PMDFs. However, besides the overall cluster, we also anal- 
yse separately the central and outer regions (according to 
the RDPs in Fig. [S}, to search for spatial variations in the 
stellar kinematics that may be related to dynamical evolu- 
tion. The deconvolved PMDFs are subsequently fitted with a 
Gaussian profile {PMDF{n) oc e-0-5((M-A)/<T)" ^j^i^j^ p^.^, 
vides the velocity dispersion (a) and the average velocity 
(ft) of the stars in the region. The field-subtracted and de- 
convolved PMDFs are shown in Figs. [S] and [S] and the rele- 
vant profile parameters are given in Table [2] The observed, 
field, and error PMDFs for the full OC sample are shown in 

App.m 

ft is clear that the deconvolved PMDFs (with dispersion 
(Jdec) are significantly narrower than the field-subtracted 
ones (aobs)- Indeed, we find adec ~ (0.1 — 0.25) aobs, but 
this ratio should be taken as an upper limit, since the ob- 
served PMDFs are not Gaussian and the fit is dominated by 
the profile core. 

The central and outer parts of NGC 1039, NGC 2516, 
and NGC 7762, have PMDFs essentially characterised by a 
single Gaussian. In some cases (NGC 2682, NGC 2477, and 
only marginally in NGC 2516) however, the observed PMDF 
has a wing towards high-PM values, which appears to be a 
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Table 2. On-sky proper motion components 



Observed Deconvolved 



Region 

(') 
(1) 


N 

(stars) 
(2) 


(mas yr 
(3) 


(mas yr -"-) 
(4) 


(mas yr""-'-) 
(5) 


(mas yr ^) 
(6) 


(kms-i) 
(7) 


(kms-i) 
(8) 


^ ohs j ^ dec 

(9) 








NGC1039 - 


— age ~ 300 Myi 


■, dQ ~ 0.5 kpc 








— 20 


204 


1.2 


> 2.4 


1.20 ±0.04 


0.45 ± 0.04 


2.8 ± 0.1 


1.0 ±0.1 


> 5.1 


5—20 


160 


a 1.2 


> 2.5 


1.23 ±0.04 


0.46 ± 0.04 


2.9 ±0.1 


1.1 ± 0.1 


^ 5.2 


0—5 


44 


1.0 


> 1-9 


1.05 ±0.04 


0.39 ±0.04 


2.4±0.1 


0.9 ±0.1 


^ 4.9 








NGC 2477 - 


— age ~ 900 Myi 


■, d© ~ 1.3 kpc 








0—20 


1030 


ai 3.3 


> 5.3 


3.44 ± 0.05 


1.12 ±0.05 


21.4 ±0.3 


6.9 ± 0.3 


^ 4.6 


0—20 


205 






9.00 ±0.50 


1.40 ± 0.40 


55.9 ±3.1 


8.7 ±2.5 




5—20 


850 


RJ 2.8 


S 4.8 


3.48 ± 0.04 


1.22 ±0.04 


21.7 ±0.3 


7.6 ± 0.3 


^ 3.9 


0—5 


232 


a 4.5 


> 7.5 


4.47 ±0.05 


1.19 ±0.05 


27.8 ±0.3 


7.4 ±0.3 


^ 6.0 


0—5 


153 






12.5 ± 0.5 


1.25 ±0.25 


77.6 ±3.1 


7.8 ± 1.6 










NGC 2516 - 


— age ~ 200 Myr, cZq 0.4 kpc 








0—20 


263 


Si 2.7 


^ 3.6 


2.72 ±0.02 


0.82 ±0.02 


4.9 ±0.1 


1.5 ± 0.1 


^4.3 


5—20 


207 


2.5 


^ 3.4 


2.56 ±0.02 


0.76 ±0.02 


4.6 ±0.1 


1.4 ±0.1 


^ 4.5 


0—5 


56 


Si 4.1 


^ 4.2 


3.89 ± 0.02 


0.99 ±0.02 


7.0 ±0.1 


1.8 ±0.1 


^ 4.2 


NGC2682 (UCAC3) — age ~ 4Gyr, ~ 0.8 kpc 


0—20 


333 


Si 2.1 


^ 3.2 


2.04 ±0.02 


0.78 ± 0.03 


7.6 ±0.1 


2.9 ± 0.1 


^4.1 


5—20 


194 


Ki 2.2 


> 3.5 


2.22 ±0.03 


0.81 ±0.03 


8.3 ±0.1 


3.0 ±0.1 


^ 4.3 


0—5 


129 


Ki 1.9 


Z 2.8 


1.77 ±0.03 


0.70 ±0.03 


6.6 ±0.1 


2.6 ±0.1 


^ 4.0 


0—5 


9 






8.20 ±0.24 


0.88 ±0.24 


30.7 ±0.9 


3.3 ±0.9 




NGC2682 (optical) 


0—15 


752 


1.4 


S 3.4 


1.75 ±0.02 


0.70 ±0.02 


6.6 ±0.1 


2.6 ± 0.1 


^ 4.9 


0—15 


113 






8.32 ±0.17 


1.13 ±0.17 


31.1 ±0.6 


4.2 ± 0.6 




5—15 


535 


1.7 


Z 4.1 


2.01 ±0.03 


0.80 ± 0.03 


7.5 ±0.1 


3.0 ±0.1 


^ 5.1 


5—15 


93 






8.69 ±0.15 


1.10 ±0.15 


32.5 ±0.6 


4.1 ± 0.6 




0—5 


217 


^ 1.3 


^ 2.4 


1.43 ± 0.02 


0.61 ±0.02 


5.4±0.1 


2.3 ±0.1 


^ 4.0 


0—5 


20 






7.86 ±0.16 


0.80 ±0.16 


29.4 ±0.6 


3.0 ±0.6 










NGC 7762 


— age ~ 2 Gyr, 


c(q ~ 1.0 kpc 








0—20 


294 


6.4 


^ 11-8 


6.70 ±0.02 


1.22 ±0.02 


30.1 ±0.1 


5.5 ± 0.1 


S 9.7 


5—20 


181 


^ 5.5 


Z 11-2 


5.69 ±0.02 


1.13 ±0.02 


25.6 ±0.1 


5.1 ± 0.1 


> 10 


0—5 


113 


Si 8.5 


^ 12.7 


8.81 ±0.02 


1.23 ±0.02 


39.6 ±0.1 


5.5 ± 0.1 


> 10 



Table Notes. Col. 2: number of member stars in region. ^ and cr derived from the fit PMDF(p) oc e 0.5((n Conversion from 

masyr"^ to kms~^ was based on the respective cluster distances (Table [TJ. 



signature of structure in the profile. This structure shows 
up in the deconvolved PMDFs as a second PM component, 
which is conspicuous especially in the central cluster region. 
The second component is characterised by a higher average 
velocity than the main one, but both have similar values of 
dispersion. 

In the central parts of NGC 2477, the second compo- 
nent corresponds to about 40% of the member stars in the 
region. This fraction drops to ~ 17% when the overall clus- 
ter is considered. For NGC 2682, this component contains 
only ~ 7% of the stars in the central region. This fraction is 
consistently the same for the UCAC3 and optical data. Ifow- 
ever, in the outer parts of NGC 2682, the second component 
is only seen in the optical data. Possible reasons for the 
difference are: the optical PMDF corresponds to a smaller 
and somewhat more interior region (5'- 15' as compared to 
5'- 20'), and contains ~ 3 times as much stars as that of 



UCAC3. The bump appears only when the observed PMDF 
is clearly non-asymmetric with respect to the average ve- 
locity, thus displaying a broad wing towards high velocities, 
e.g. NGC 2477 and NGC 2682. On the other hand, there is 
no bump emerging from the essentially gaussian (observed) 
PMDFs of NGC 1039, NGC 2516, and NGC 7762 (Figs. E] 
and [6]). In addition, it is interesting to note that the bump 
in the central region of NGC 2682 appears almost identically 
in PMDFs built with independent data sets. 

Another issue is to what degree incompleteness in 
crowded regions - and the more difficult measurement of 
PM components for faint stars - affect the PMDFs. Addi- 
tionally, could the bump be related to incompleteness? We 
use NGC 2477, the most distant and populous OC of our 
sample (thus, the most prone to suffering from incomplete- 
ness) to examine this point (App. IB1[) . Since our analysis 
depends on the number of member stars (especially at the 
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NGC2477 



NGC2516 



NGC7762 




5 10 15 20 25 5 10 _|15 10 20 30 

H„, (mas yf ) 

Figure 6. Same as Fig. [5] for the remaining clusters. 



central region), we selected J — 12.5 as the boundary be- 
tween bright and faint stars. At the distance of NGC 2477, 
this boundary corresponds to a stellar mass of m ~ 1.6 Mq. 
The bright and faint deconvolved PMDFs built for the cen- 
tral region are similar (Fig. IB3|l , consistently presenting the 
high- velocity bump at /7 ~ 12 mas yr~^. The only signifi- 
cant difference is that the faint (less massive) star PMDF is 
shifted ~ l.Smas yr"'^ towards high values of jl with re- 
spect to that of the bright stars. A similar shift occurs 
for the bright and faint star PMDFs in the outer region 
(without the bump). Although relatively small, the shift 
Afi ~ l.Smas yr~^ ~ 8kms~^ between the bright and faint 
PMDFs might suggest sUghtly different kinematics for stars 
in different mass ranges. This experiment also suggests that 
incompleteness and PM measurements of faint stars - at 
least to the level available in UCAC3 - are not critical for 
the PMDFs. 

Additionally, one might ask whether the bump may 
come from residual, i.e. unaccounted for field contamination. 
It is true that, given the statistical way we decontaminate 
the clusters (Sect. (3)), some field-star contribution might per- 
sist in the subtracted PMDFs. However, both NGC 2477 and 
NGC 2782 are located in the third Galactic quadrant and at 
high Galactic latitudes, which by itself minimises contami- 
nation (as can also be seen in Figs. lBl1 and lB2] ). Thus, any 
residual contamination should be minimum, which would 
contradict the fraction of stars composing the bump, ~ 40% 
in the central region of NGC 2477 and « 7% in NGC 2682. 
The above arguments suggest that the bump is related to 



a cluster kinematic property, but we cannot definitely rule 
out the possibility that it might be an artefact of the RL 
deconvolution, and/or some residual contamination by stars 
with peculiar PM components. 

There is no direct interpretation for the additional 
bump seen in the deconvolved PMDFs of NGC 2477 and 
NGC 2682. Assuming that it is physical, one possibility is 
that the double peak may arise from a merger of two clus- 
ters (as s uggested by the referee, Thi j s Kou wenhoven). In 
this sense, Ide Oliveira. Bica fc Dottoril (|2000l ) carried out N- 
body simulations of cluster encounters, studying long-term 
structural changes up to ~ 1 Gyr. They found that the clus- 
ters may coalesce at such ages, but until then, the presence 
of the two clusters can still be traced by means isophotal 
distortions and elliptici ty variations, as observ ed in model 
and actual clusters (e.g. Ide Oliveira et al.|[200 ^). In this sce- 
nario, internal differences in kinematics might persist too, 
producing different signatures in the deconvolved PMDFs. 

Alternatively, the high-average velocity component may 
be related to mass segregation, in which a fraction of the 
stars collectively migrate along the radial direction (over a 
time-scale of a relaxation time) within a star cluster. In the 
cases dealt with here, it occurs only in the two most pop- 
ulous DCs, NGC 2477 and NGC 2682. Maybe it cannot be 
detected (by the present approach) in the other OCs because 
they are less populated (i.e., possibly the same reason why it 
is seen in the outer parts of NGC 2682 with the optical data, 
but not with UCAC3). Fi nally, working wi t h hist ograms and 
a different PM-data set, iBica fc Bonattd l|2005h raised the 
possibility that the high-velocity component in NGC 2682 
might be relat ed to the presence of binar i es. Ho wever, given 
the findings of iKouwenhoven fc de Griid (|2008l '). this possi- 
bility seems the least probable. In any case, a definitive solu- 
tion for the bump nature would require detailed simulations 
of the internal cluster dynamics (including cluster merger, 
mass segregation, and binarity), a task that is beyond the 
scope of the present paper. 

The velocity dispersions derived from the deconvolved 
PMDFs (Table ^ of NGC 1039 and NGC 2516 (and, to a 
lesser degree, NGC 2682), are consistent with those of (ap- 
proximately virialized) OCs of ~ 10'' Mq. On the other 
hand, those of NGC 2516 and NGC 2477 appear to be exces- 
sively high for OCs of a similar mass scale. However, such 
large dispersion values can be partly explained by an obser- 
vational limitation related to distance, since both OCs are 
the most distant of the sample. For a limited observational 
time-base, the PM determination for a distant OC will pri- 
marily detect the high-PM components, thus leading to a 
broad PM profile and high values of velocity dispersion. Un- 
der similar conditions, a nearby OC, on the other hand, will 
also have part of the low-PM components detected, thus im- 
plying a lower velocity dispersion. This effect is present in 
our analysis, as can be seen by the correlation between CTos 
and distance from the Sun (Fig. [7] - left panels) . An addi- 
tional consequence of the distance-related effect would be a 
shift towards high values of the on-sky average velocity (fi) . 
Again, this shift would increase with the distance from the 
Sun, and this relation is also present in Fig. [T] Despite this 
effect, the average velocity of the bump stars in NGC 2477 
(/i ~ 78kms~^) is exceedingly high for a ~ 10^ M© clus- 
ter, to the point that these stars - if they really belong 
to NGC 2477 - are not gravitationally bound to the clus- 
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Figure 7. Dependence of the deconvolved profile parameters 
on distance from the Sun (left panels) and cluster age (right). 
The central (empty circles) and outer (filled squares) regions are 
shown. Given the rather limited range of values oi aabs/'^d£c (Ta- 
ble [2}, similar relations hold as well for the observed values of p, 
and <Jos- 



ter. This suggests that the high-velocity bump may be an 
artefact produced by the RL deconvolution (when applied 
to asymmetric profiles) or, least probably, unaccounted for 
field contamination. 

Given the above caveat, it is not possible to disentangle 
a physical relation between aos and cluster age (Fig. [7]- right 
panels) from the observational limitation. 



phology for establishing the colour and magnitude ranges of 
the probable member stars, (ii) these stars are used to build 
the RDP, which provides the cluster structural parameters 
and allows to define the comparison field, (Hi) we define a 
grid of PM bins of variable size that spans the full range of 
values of /ia cos((5), fis, and /ioa, (iv) considering the cluster 
and field stars separately we compute the probability that 
the PM measurements of a given star corresponds to any 
bin, (v) the field-star PMDF is subtracted from the cluster 
PMDF and, (vi) we take the intrinsic PM-error distribu- 
tion function as the PSF to be used in the Richardson-Lucy 
deconvolution approach. 

The main result of our approach is that the deconvolved 
PMDFs are well represented by Gaussians with dispersions 
lower than the observed ones by a factor of 4-10. Besides 
the main component, the deconvolution revealed structure 
in the profiles of NGG 2477 and NGC 2682 in the form of 
a second - and less populous - distribution shifted towards 
higher average velocities, which may originate from cluster 
merger, large-scale mass segregation or, least probably, bi- 
naries. The secondary bump in NGG 2477 consistently ap- 
pears in PMDFs built with stars brighter and fainter than 
J = 12.5, which suggests that it is not related to incom- 
pleteness and/or faint star-PM measurement. However, we 
cannot exclude the possibility that the bump is an arte- 
fact produced by the RL deconvolution when applied to 
strongly asymmetric profiles such as those of NGG 2477 and 
NGC 2682. NGC 1039 and NGC 2516, the nearest OGs anaL 
ysed here, end up with deconvolved dispersions compatible 
with those expected of bound OGs of ~ 10"^ Mq . We also 
detect an increase of the velocity dispersion and average ve- 
locity with distance from the Sun, which is probably due to 
a similarly limited time-base used for measuring the proper 
motions among different star clusters. 

In recent years our group has given particular attention 
to the investigation of OCs by means of analytical tools that 
produce field-star decontaminated GMDs and RDPs. These 
tools have proved essential for a constrained analysis of OGs 
characterised by a range of parameters (age, distance from 
the Sun, reddenin g, etc), and located on a wide variety of 
environments (e.g. Bonatto fc Bical l2010bl : iBonatto &: Bical 
l2007l : iBonatto et al.ll2006l ). The present paper links the clas- 
sical GMD and RDP analyses with a novel approach for 
dealing with proper motions and the respective errors. 



6 SUMMARY AND CONCLUSIONS 



A crucial point in understanding a cluster's dynamical stage 
is the derivation of kinematical parameters for its member 
stars, usually by means of proper motions obtained in public 
databases. However, uncertainties associated with ground- 
based proper motion measurements are usually large, and 
their effect should be properly taken into account when 
building PM profiles for determining the velocity dispersion. 

In this paper we investigate the above issue us- 
ing the relatively nearby and populous open clusters 
NGC 1039 (M34), NGC 2477, NGC 2516, NGG 2682 (M67), 
and NGC 7762 as test cases. Their PM components have 
been obtained in UCAC3. 

Rather than working with PM histograms, we build 
PMDFs for the cluster and field stars, taking the la-PM un- 
certainties into account. In short, (i) we use the GMD mor- 



ACKNOWLEDGEMENTS 

We thank the referee, Thijs Kouwenhoven, for interesting 
comments and suggestions. We acknowledge support from 
the Brazilian Institution GNPq. This publication makes use 
of data products from the Two Micron All Sky Survey, which 
is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Centre/California In- 
stitute of Technology, funded by the National Aeronautics 
and Space Administration and the National Science Foun- 
dation. We also employed the WEBDA database, operated 
at the Institute for Astronomy of the University of Vienna. 



Deconvolved OC-velocity dispersion 9 



REFERENCES 

Bi H. & Boerner G. 1994, A&AS, 108, 409 
Bica E. & Bonatto C. 2005, A&A, 431, 943 
Bonatto C, Bica E. & Girardi L. 2004, A&A, 415, 571 
Bonatto C., Bica E., Ortolani S. & Barbuy B. 2006, A&A, 
453, 121 

Bonatto C. & Bica E. 2007, MNRAS, 377, 1301 

Bonatto C. & Bica E. 2009b, MNRAS, 394, 2127 

Bonatto C. & Bica E. 2010a, A&A, 516, 81 

Bonatto C. & Bica E. 2010b, A&A, 521A, 74 

Chen L., de Grijs R. & Zhao J.L. 2007, AJ, 134, 1368 

Girardi L., Bertelli G., Bressan A., Chiosi C., Groenewegen 

M.A.T., Marigo R, Salasnich B. & Weiss A. 2002, A&A, 

391, 195 

Goodwin S.P. & Bastian N. 2006, MNRAS, 373, 752 

de Grijs R., Goodwin S.P., Kouwenhoven M.B.N. & 

Kroupa, P. 2008, A&A, 492, 685 
Kouwenhoven M.B.N. & de Grijs R. 2008, A&A, 480, 103 
Lada C.J. & Lada E.A. 2003, ARA&A, 41, 57 
Lamers H.J.G.L.M., Gieles M., Bastian N., Baumgardt H., 

Kharchenko N.V. & Portegies Zwart S. 2005, A&A, 441, 

117 

Lucy L. 1974, AJ, 79, 745 

de Oliveira M.R., Bica E. & Dottori H. 2000, MNRAS, 311, 
589 

de Ohveira M.R., Dutra CM., Bica E. & Dottori H. 2000, 

A&AS, 146, 57 
Pavani D.B. & Bica E. 2007, A&A, 468, 139 
Piskunov A.E., Schilbach E., Kharchenko N.V., Roser S. & 

Scholz R.-D. 2007, A&A, 468, 151 
Richardson W.H. 1972, J.Opt.Soc.Am., 62, 55 
Richer H.B., Dotter A., Hurley J., Anderson J., King I., 

Davis S., Fahlman G.G., Hansen B.M.S. et al. 2008, AJ, 

135, 2141 

Skrutskie M.F., Cutri R., Stiening R., Weinberg M.D., 
Schneider S.E., Carpenter J.M., Beichman C, Capps R. 
et al. 2006, AJ, 131, 1163 

Spitzer L. 1987, in Dynamical Evolution of Globular Clus- 
ters, Princeton, NJ, Princeton University Press, p. 191 

Vio R., Bardsley J. & Wamsteker W. 2005, A&A, 442, 397 

Yadav R.K.S., Bedin L.R., Piotto G., Anderson J., Cassisi 
S., ViUanova S., Platais I., Pasquini L. et al. 2008, A&A, 
484, 609 

Zacharias N., Monet D.G., Levine S.E., Urban S.E., Gaume 

R. & Wycoff G.L. 2004, AAS, 205, 4815 
Zacharias N., Finch C, Girard T., Hambly N., Wycoff G., 

Zacharias M.I. , Castillo D., Corbin T. et al. 2010, AJ, 

139, 2184 



APPENDIX A: SOME ASPECTS RELATED TO 
THE RL DECONVOLUTION 

Al PMDF resolution 

For preserving the profile resolution especially around the 
PMDF peak (where the number density of stars is high) 
and, at the same time, keep acceptable error bars (especially 
towards the PMDF wings where the number density is low), 
we adopt a variable distribution of bin sizes Spm , where pm 
represents any component among {fiaCos{S), ns, ^os) for 
building the PMDFs. 




1 10 100 1 10 100 



iterations iterations 

Figure Al. Top: variation with the number of deconvolution 
iterations of the on-sky velocity dispersion Oos (left) and average 
velocity p, (right). Bottom: fractional variation of aos (left) and 
X (right). The curves correspond to the 0' — 20' PMDFs, except 
for the optical PMDF of NGC 2682 buih for 0' - 15'. 

After determining where the PMDF peak occurs, we 
use (5pm = 0.25 mas yr~^ for shifts with respect to the peak 
in the range |Apm| = — lOmas yr~^, Spm ~ Imas yr~^ 
for |Apm| = 10 — 30masyr~^, Spm = 5masyr^^ for 
Ap„i| = 30 — 50masyr~^, and 5pm = lOmas yr~^ for 
Apm > 50 mas yr~^ and Apm < —50 mas yr~^. Addition- 
ally, this procedure also has the advantage of reducing com- 
putation time. 



A2 Deconvolution convergence 

The rather slow convergence rate of the RL deconvolu- 
tion method has been extensively discussed in the literatur e 
Ce.g. lBi & Boerneilll994l : IVio. Bardslev & Wamste"kei1l2005l '). 

However, it is easy to implement, reliable, and the quality 
of the outputs can be followed at each iteration. 

We follow in Fig. lAll the changes on aos and /i with the 
number of deconvolution iterations {Nu). While the average 
velocity (panel b) changes little over the 200 iterations ap- 
plied here, o-os, on the other hand, decreases systematically 
with Nit (panel a) but with a rate that begins to flatten 
for Nit ^ 100. Clearly, aos could decrease somewhat more 
for Nit > 200. However, as implied by the rate of change 
faoa ~ Aaos/o'os that, for Na = 200 has decreased to 
faoa lS 0.003 (panel c), it would take several hundred more 
iterations to produce a significant change in aos- Finally, 
the chi-square (Sect.|4]) fractional variation — Ax/x ^Iso 
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Figure Bl. The panels show the observed, field, and intrinsic 
error (PSF) on-sky PMDFs built for different cluster regions: cen- 
tral (top), outer (middle) and overall (bottom). The number of 
member stars in each region is indicated. 



presents a steady decrease with Nu (panel d), dropping to 
~ 0.003 (or a ~ 0.3% variation) for Nu = 200. 
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Figure B2. Same as Fig. [Bl] for the remaining OCs. 

the "full-magnitude" range PMDF, both for the central (left 
panel) and outer (right) cluster region. 



APPENDIX B: OBSERVED AND FIELD PMDFS 

We show in Figs. IBll and IB2I the observed on-sky PMDFs 
for all clusters of our sample, built for the central, outer, 
and overall cluster regions. Also shown are the PMDFs cor- 
responding to the intrinsic error distribution (PSF) and the 
field. It is interesting to note the relative contribution of the 
field stars PM among the different regions of the same clus- 
ter and among the full set of clusters. This shows that the 
field contribution must be taken into account in the analy- 
sis. Also, and perhaps more importantly is the PSF width 
that, in some cases, is almost as broad as the cluster PMDF. 



Bl PMDFs in different magnitude ranges 

In Fig. lB3l we examine the incompleteness/faint stars issue in 
NGC 2477, which should be important especially for the cen- 
tral part of the cluster. Starting with the CMD of NGC 2477 
(Fig. [2|, we build PMDFs separately for stars brighter and 
fainter than J = 12.5. This magnitude boundary is adequate 
for characterising different types of stars and, at the same 
time, keeping a reasonable number of stars in each magni- 
tude range. The bright and faint PMDFs are compared to 
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Figure B3. Deconvolved PMDFs for the central and outer re- 
gions of NGC 2477, separated by magnitude ranges. The bright 
(solid line) and faint (dashed) PMDFs are compared to the "full- 
magnitude" (shaded) PMDF. The number of stars composing the 
PMDFs is indicated. 



